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Asymptotic behaviour of the probability density in one dimension 
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We demonstrate that the probability density of a quantum state moving freely in one dimension 
may decay faster than t~ l . t~ 2 and t~ 3 dependences are illustrated with analytically solvable 
examples. Decays faster than t -1 allow the existence of dwell times and delay times. 
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T D (a,b) = dt dx\i>{x,t)\ 2 , (1) 

J —oo J a 

and delay times for one dimensional motion, the later being defined as the difference between the dwell time with and 
without an interaction potential. El n 
The classical counterpart of ([!]) isO 



7— I \ The study of the asymptotic, long-time behaviour of wave packets, and related quantities such as decay and 
■ survival probabilities has a long history. Recently UnnikrishnanHH has proposed several examples of decay of the 
wavefunction which are slower than the £ -1 / 2 decay that follows from straightforward dimensional analysis. In 
particular, Unnikrishnan has put forward behaviours of the type t~ x / 2+ll ,Q < [i < 1/4 and t~ 1 / 2 \at. In this paper 
Oh' we complement this study by giving examples of decays which are faster than those considered in Refsjl],^. We give a 
procedure for obtaining systematically different kinds of faster decays, which we illustrate with two explicit examples, 
and then discuss the failure of the arguments that suggest that only the < -1 / 2 (or slower) decay canlp,ke place. 
An important consequence of the existence of faster decays is the possibility to define dwell timesjtl 
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O" T D (a,b) classical = I dt I dxg(x,t), (2) 

Qh where g(x, t) is the probability density of an ensemble of independent particles. rrj(a, b) c i aS sicai is the average over the 
ensemble of the time that each particle trajectory spends between a and b. In other words, this is an average "dwell" 
or "sojourn" time in the selected region. Even though the interpretation of ([j]) is not as straightforward, since the 
individual member of the quantum ensemble is not associated with a trajectory in the standard interpretation, the 
quantum dwell time is in any case recognized as an important characteristic quantity of the wave packet J3 Similarly, 
the delay time is, together with the transmittance, one of the basic quantities that summarize the essential aspects 
of wave packet scattering, also in the stationary limit, to can be written in several ways, in particular as 

T D (a,b)= f dtP ab (t) = (iP(t = 0)\T D \i(;(t = 0)), (3) 

_ _ I J — oo 

where P a b{t) — f dx ^(x, i)| 2 , Tjj is the sojourn time operator, 

/oo ^ ^ 

dte lHt / h D(a,b)e- lHt/h , (4) 
-oo 

H is the Hamiltonian, and D(a, b) is the projector onto the selected space region, 

D(a,b) = [ dx\x)(x\. (5) 



Unless P a b(t) decays faster than t , the dwell time will diverge. The existence of a potential function leads generically 
to an asymptotic decay ~ t~ 3 , as discussed in Ref. || (see also Ref. || and references therein). The situation is however 
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rather different for the case of free motion, with Hamiltonian Hq, for which the sojourn time operator (Q) takes in 
momentum representation the form 

/OO T 
dp^\p)(p\D\ap)(ap\. (6) 
,1-1 -°° >P> 



Note the presence of a factor that will make the dwell time divergent unless the state vanishes at p = 0. The 

fact that generically the dwell time for free motion will diverge also follows from the dependence on i 1 / 2 of the 
corresponding propagator, 

(x\e- iSot / H \x') = (™\ 1/2 jrn{*-* 072^ ( 7 ) 

\ihtJ 

since this leads to a i -1 dependence of the probability density \ip(x,t)\ 2 in a generic case. 

There exists an equivalent scaling argument to understand this fact: in the generic free motion case there are no 
dimensional magnitudes other than fi and m present, and therefore no relevant magnitude with dimensions of time 
can exist. Since P a b{t) is, dimensionally speaking, an inverse time, it has to scale with i -1 in the absence of such 
a quantity. In the case of interaction, the interaction length l Q comes into play, thus modifying this generic scaling, 
since the magnitude ml^/fi is a characteristic interaction time. 

On the other hand, if the initial data is such that a characteristic length or a characteristic momentum exists, and 
it cannot be eliminated by translations or Galilean boosts, it is likely that the large time asymptotics will be modified. 
An alternative way of looking into this possibility is by noticing that, using Eq. (0), an integral over x' is required 
to obtain ip(x,t), which entails that the asymptotic i -1 / 2 dependence of the wavefunctions is not guaranteed in all 
cases: the limit ipi-> oo and the integral do not necessarily commute. In particular, Unnikrishnan has found examples 
of slower decay,™ and we shall justify and illustrate decays faster than i -1 when the momentum amplitude (p\ip) 
vanishes at p = 0. |— ■ 

The large time asymptotics of i(j(x, t) will depend on the behavior near the "critical point"l!3 p = 0, 

/■oo 

#M) = W(*)>= / dpix^e-^/^ipim), (8) 



whereas larger momenta are responsible of the short time behaviour, as happens when considering a classical ensemble 
too. To determine the large time behavior we shall assume that (p\ip(0)} may be expanded around the origin (this 
is not the case for the "Cauchy state" chosen by Unnikrishnan to show slow decayBEl - in fact, this state presents a 
logaritmic singularity at the origin of momenta). Let us assume in particular that the state belongs to the subspace 
of positive momenta, and that it vanishes at p = 0. Without loss of generality we set x = 0, 

/>OC 

i/j(0,t) ~ h- 1 ' 2 I dpe~ lp2t / 2mh [ap + bp 2 + ...}. (9) 

Now make the change of variables 
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u=p/f, /=(l-t)V(mV*)- (10) 
On deforming the contour to the steepest descent path in the fourth quadrant (where u is real), we obtain 

(x = 0|VW) ~ h- 1 ' 2 ! / due~ u \afu + b{fu) 2 + ...], (11) 
Jo 

which leads to a large time dependence \tp(x = 0, t)\^-f= 0[t~ 2 ] for a ^ 0, or 0[i -3 ] for a = 0, b ^ 0. 
An example of the later case is given by the stateO 

0#(0)) = C(l - e -"f 2 /R 2 ) e -« 2 (p-Po) 2 /fi 2 -^o/n ( p ) ; ( 12 ) 

for which a = 0, and 6^0. (C is the normalization constant and is the Heaviside "function" .) The time dependence 
may be obtained explicitly from (||) , 

WW) = 47rW3 * I Zv^ (XT^T 2 — /• (13) 




FIG. 1: dln\(x\ip(t))\ 2 /dlnt versus dint for two different wave packets: one of them is given by (^2(), (solid line), and the other 
one is a a Gaussian wave packet, <p|V(0)} = C'e~ s (p - po) /n -w^o/* (dashed line), where C is the normalization constant. 
The parameters are po = 1> xo = —10, a = 0.5, 5=1, x — 0, m — 1 (all quantities in atomic units). Note the asymptotic 
dependences of the probability densities: t -3 and t -1 respectively. 



Here w(z) = e~ z erfc(— iz), fco = pa/fr, A = S 2 +iht/ (2m), and g = i(x — xq) + 2ko5 2 . From the asymptotic properties 
of the w function, see the appendix, it follows that \ip(x,t)\ 2 — 0[t~ 3 ] as illustrated in Figure 1, where we have also 
depicted the t~ 1 decay corresponding to a minimum-uncertainty product state. 

The intermediate decay t~ 2 in the probability density is obtained for example for the initial wavefunction 

N 

V>0,o) = - — —2 , 

(x + la) 



where a is real and positive, and N = ^2a 3 /TT is the normalization factor. As a matter of fact, this is a direct 
counterexample to the statement of Unnikrishnan that spatial decay faster than |rr:| 1 necessarily entails temporal 
decay of the form t~ x l 2 for the wavefunction. The Fourier transform of this initial wave packet is 

/a\3/2 



ipip, 0) = -2 f|) <d(p)pe- a P/ h 



from which it is straightforward to conclude the decaying behaviour mentioned above. The explicit expression for all 
times can be computed in terms of 1-F2 hypergeometric functions. 

Another example with t~ x decay in the wavefunction is given by considering the following initial wavefunction in 
momentum representation: 

^feo) = - T7 ^ 7I e(^ e - 2 ^. 

The solution of Schrodinger's equation with this initial condition, now written in coordinate space, is 

7T J '*(«pt + rah) - 



where 
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Again using the series expansion of the w function, one can easily check the asymptotic behaviour mentioned above. 
In fact, in this particular example it is easy to derive an explicit expression for din | (x\ip(0, t))| 2 /c?lni, which results 
in -2/(1 + m 2 h 2 /p 2 t 2 ). 

In conclusion, we have shown that-|tlie asymptotic, large time decay laws of the free-motion one-dimensional wave- 
function described by Unnikrishnantra (as i -1 / 2 or slower) are not exhaustive. Faster decays (i -1 and £~ 3 / 2 ) have 
been demonstrated with analytically solvable examples and may be expected for wave packets that vanish at zero 
momentum. These faster decays lead to finite dwell times and allow to give a meaning to the "lifetime" (delay timpl 
matrix of Smithtj in one dimension in terms of the difference of dwell times with and without interaction potential.Ej 
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APPENDIX A: PROPERTIES OF m-FUNCTIONS 



The w— function is an entire function defined in terms of the complementary error function a w(z) 
e~ z exfc(— iz). Its integral expression is 



1 f e~ u 

w{z) = — / du, (Al) 

ITT r U — Z 



where T_ goes from — oo to oo passing below the pole at z. From (Al) two important properties are deduced 



w{-z) = 2e~ z -w{z) (A2) 

and 

w(z*) = [w(-z)}*. (A3) 
The derivative of the w-function can be expressed in terms of the w-function itself, 

2i 

w'(z) = -2zw(z) + —=. (A4) 

Finally, w(z) has the series expansion 

-W'Sflfl)- < A5 » 
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